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Abstract 

Bondi type hydrodynamic accretion of the surrounding matter onto Kerr black hole with an arbitrary rotational 
parameter is considered. The effects of viscosity, thermal conductivity and interaction with radiation field are 
neglected. The black hole is supposed to be at rest with respect to matter at infinity. The flow is adiabatic 
and has no angular momentum. The fact that usually in astrophysics substance far from the black hole has 
nonrelativistic temperature introduces small parameter to the problem and allows to search for the solution as a 
perturbation to the accretion of a cold, that is dust-like, matter. However, far from the black hole on the scales 
of order of the radius of the sonic surface the expansion must be performed with respect to Bondi spherically 
symmetrical solution for the accretion on a Newtonian gravitating centre. The equations thus obtained are 
solved analytically. The conditions of the regularity of the solution at the sonic surface and at infinity allow 
to specify unique solution, to find the shape of the sonic surface and to determine the corrections to Bondi 
accretion rate. 
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1 Introduction 

The problem of the matter accretion onto gravitating objects is the classical problem of astrophysics (Zeldovich 
& Novikov 1967; Shapiro & Teukolsky 1983). The accretion is widely engaged for the explanation of various 
kinds of astrophysical phenomena. It is believed that this process is basic to understanding the nature of 
active galactic nuclei, quasars, galactic X-ray sources, the formation of extragalactic jets and jets from some 
galactic objects of stellar masses (i.e., Begelman, Blandford & Rees 1984; Camenzind 1990; Lipunov 1992). 
Usually, when considering accretion with astrophysical implications, the authors either worked in the frame 
of a disk accretion model or considered spherically symmetrical case, which enabled them without excessive 
details to estimate the accretion rate and the emerging radiation due to the heating of infalling matter (Shapiro 
& Teukolsky 1983). Bondi (1952) was the first who considered hydrodynamic isoentropic stationary flow for 
polytropic equation of state and showed that, being subsonic at infinity, flow must cross the sonic surface and 
become supersonic near the black hole, the accretion rate being determined by the requirement of smoothness 
of the solution at the sonic surface. 

Present paper investigates the adiabatic stationary accretion onto a rotating black hole. Let us first note that 
the luminosity of the infalling matter itself is much less than the Eddington limit in quasispherical case (Shapiro 
& Teukolsky 1983), so, despite that our approach does not involve interaction with photons, it can be applied for 
the accretion by an isolated black hole embedded in interstellar gas or molecular cloud. Probably Shapiro (1974) 
was the first who considered accretion onto a Kerr black hole. He noticed that in the case of zero temperature the 
situation is similar to the accretion of dust, so the particles would follow geodesic lines in Kerr space-time and 
flow lines would simply coincide with them. As it can be verified directly from the equations for the test particle 
motion in Kerr geometry (e.g., Novikov & Frolov, 1986), a particle, resting with respect to the black hole at 
infinity, will move purely radially (in Boyer-Lindquist coordinates) and experience dragging in (j) direction with 
Lense-Thirring velocity, while any additional impact to it in radial direction causes it to deflect in 6 direction. 
Therefore, we already have the accretion pattern for cold matter. One should also expect that the slight heating 
of the gas far from the black hole, such that the gas molecules thermal velocity remains to be much smaller than 
the velocity of light, would not alter the flow picture drastically. Having had noticed all that, Shapiro came 
to calculation of the emerging bremsshtralung radiation from almost cold infalling material. Further step was 
undertaken by Petrich, Shapiro & Teukolsky in 1988. They found analytical solution for arbitrary fast rotating 
Kerr black hole moving through the surrounding gas along its rotational axis, but for equation of state p = Q 
with the sound velocity equal to the light velocity everywhere, which looks very artificial. 

Recently Beskin & Pidoprigora (1995) (further referred to as BP) considered the accretion for polytropic 
equation of state by slowly rotating black hole. They used small ratio of the rotational parameter a to the black 
hole mass M. as a small parameter to linearize hydrodynamic equations and to search for the solution as a small 
perturbation to the spherically symmetrical accretion by Schwarzschild black hole. The aim of the present work 
is to find the solution for the Kerr black hole with an arbitrary value of a. As in BP the black hole is supposed 
to be at rest with respect to the gas at infinity, flow is adiabatic and without angular momentum with respect 
to the black hole at infinity, matter has uniform temperature and density far from the black hole. Under these 
conditions flow must be axisymmetrical with the symmetry axis being the rotational axis of the black hole. 

For solving the problem formulated we use the equation of flow lines equilibrium (stream equation), which 
is nonlinear differential equation of mixed type in flux function $ and is of Grad-Shafranov type. Its general 
form for Kerr metric has been derived in Beskin & Pariev (1993). In section 2 basic quantities describing the 
flow are introduced and the equilibrium equation is formulated. As it is pointed above this equation must have 
solution corresponding to the accretion of the cold matter having zero pressure {p = 0). It occurred to be 
really true as have been verified in BP. The fact that usually in astrophysics substance far from the black hole 
has nonrclativistic temperature introduces small parameter to the problem — the ratio of the temperature at 
infinity to the rest mass of a particle, and allows us to search for the solution as a perturbation to the accretion 
of a cold matter. Stream equation can be linearized with respect to this parameter and the linearized equation 
occurred to be solvable analytically. This has been done in section 3. However, the solution for cold matter 
cannot be considered as a zero order solution for the expansion on the distances from the black hole comparable 
to the radius of the sonic surface (which is much greater than the radius of the event horizon), because the 
appearance of the sonic surface itself is essentially connected with the nonzero temperature of the matter. In 
section 4 we have constructed solution of the expanded stream equation in that region using as a zero order 
approximation the solution for Bondi accretion flow onto Newtonian gravitating centre of mass M.. Matching 
of the two solutions allows to determine the accretion pattern uniquely. 

Construction of the analytical shock-free solution for the hydrodynamic accretion onto Kerr black hole leads 
us to the generalization for rapid rotation of the result of the work BP that neither the rotation of a black hole 
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nor its slow motion change drastically the nature of the accretion from the spherically symmetrical smooth flow. 

2 Basic equations 

Let us consider an axisymmetrical steady-state flow in the Kerr space-time with the metric in Boyer-Lindquist 
coordinates x° =t,x^ = r, = 9, = 0: 



where 



= = 0, P^ = -u: = (3) 



ds^ = -a^df + gikidx' + (}'dt){dx^ + p'^dt), (f ) 

a = |^/A (2) 
is the lapse function (a = A = at the event horizon), 

2aMr 

LD is the angular velocity of Lense-Thirring precession, for all i 7^ fc gik = 0, 

grr = P^/A, gee = p'^, g^^ = w"^ , (4) 

and 

A = + - 2Mr, = r'^ + cos^ 6, 

= {r^ + a^f -a^/:\s\T?e, w^-sinO. (5) 

P 

Here M is the mass of a black hole, a = J /M is the specific angular momentum of a black hole. Throughout 
all the paper we will adopt the notations used in Beskin & Pariev (1993) and in BP. We will use units where 
c = 1 and G = 1, that is all physical quantities have dimensions of some power of the unit of length. We follow 
the '3+1'- splitting of the space-time described in the book by Thorne et al. (1986). All physical quantities 
considered here are measured by Zero Angular Momentum Observers (ZAMO), rotating around the black hole 
at a constant radius r with Lense-Thirring angular velocity d(f>/dt = uj. All 3-dimcnsional vectors and tensors 
are referred to 3-dimensional 'absolute' space with metric gik @) orthogonal to the world lines of ZAMOs. It is 
convenient to use the following orthonormal triad in this 3-space 



Ad Id Id 

Axial symmetry and time independence allow us to introduce flux function $(r, 6) defined as 

"""^^ = ^(^* e^)- (6) 

Here n is the total particle density in the comoving reference frame, up is the poloidal (that is in the directions 
of eg and ) component of the four velocity u^. Due to the definition of $(r, 0) by formula (||) the continuity 
equation V-(anu) = is satisfied automatically. Because of (uV$) — the flux function ^{r,0) is constant 

along each stream line, and the quantity j d^^ — j an urfS is the total flux of the particles through an area S. 

s s 

To avoid singularities on the rotational axis 9 = Q flux function must satisfy two conditions *&|g^g — and 



9$ 



'do 



= 0. Note also that the value — ^Ig^^ is the total accretion rate (with positive sign). 



Particularly, in the case of spherical symmetry $ = *i>o(l ~ cosO), and the accretion rate is — 2$o. 
Following the tradition we adopt polytropic equation of state 

p^k{s)n^, (7) 

where p is the gas pressure, 1 < F < 5/3 is the polytropic index, and fc(s) is the quantity depending only on 
the entropy per particle s. The flow is adiabatic, therefore, s is constant along stream lines $(r, 9) = const, so 
s = s(<&). Then, entalpy per particle p = (p + gm)/n (here is the specific internal energy of the gas) is 

p^m+^^k{s)n^~\ (8) 
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where m is the mean rest mass of a particle. Accordingly, the sound velocity is 



dp 
dn 



rfc(s): 



(9) 



All thermodynamic quantities can be expressed via the functions in s and n only. Moreover, in our problem gas 
temperature and density are uniform on large distances from the black hole, so the entropy s is constant along 
all stream lines, that is constant everywhere s = const. Therefore, all thermodynamic quantities effectively 
depend on only one variable and can be expressed via each other given the value of s. 

The equations of motion are provided by four components of the energy-momentum conservation low 



T 



0. 



(10) 



The 0-component of (|T^) leads to the energy conservation along each stream line (Bernoulli integral), which is 
manifested by that the following quantity E depends only on $ 



(11) 



1+ M A + U] 



is the Lorentz factor of the flow. The 0-component of ( p^ leads to the conservation of 



Here 7^ 

the z-component of the angular momentum along each stream line 

L($) = nwu,. 



(12) 



For the problem in hand £($) = 0, which means = 0. Hence, the gas corotates with ZAMOs with the 
angular velocity just equal to lu as seen by an infinitly distant observer. At infinity a = 1, 7 = 700 — 1, and 
M = Moc = const, therefore, one can see from ( pT| ) that E{^) — E — const along all stream lines and, hence, 
everywhere. Now we can express all thermodynamic quantities and the velocity u by means of <i>(r, 9) and 
constants E and s only. Using definition (|^) and obvious relation 7^ = 1 + Up and bearing in mind that L = 
one can obtain the following equation 



m^E^ = /i^ Asin^ei 



(13) 



E, 



which, being combined with the formula (^) for /i, provides an implicit expression for n in terms of |V<i>p 
and s. Then, all other characteristics of the flow can be obtained from (||) and (^). 

The remaining two space poloidal components of the equation ( p^ ) should be decomposed into parts along 
the stream line, i.e. perpendicular to V*!", and perpendicular to the stream line, i.e. parallel to V$. Component 
along the stream line leads to the adiabatic condition s — s($) and gives no new information, because even the 
form of the T^^ for ideal fluid itself involves the conservation of the entropy (viscosity and thermoconductivity, 
if any, are taken into account by including additional terms in T^^\ Component perpendicular to the stream 
lines results in the equation of stream lines equilibrium — stream equation (or the Grad-Shafranov equation). 
With the terms containing second derivatives written out explicitly it looks as follows 



A 



92$ 



d 



smf 



2A 



9$ a$ .92$ 



1 9$ 



de V sin 61 86 



1 



A(9$/5r)2 + {d<^/dey 



A2 ^ 



dr 



32$ 

Q^2 



dr 86 drde ' \de 
1 1 V 2Tr^n^ p2 



92$ 
9^ 



1 OA a$ 

2 dr dr 



'd6 



ld_ 

2dr 



A2 a$ 

p2 dr 



(14) 



where 



2p2 d6 de 



D = -1 



|V$|2 ^2 



V$ • [E^V{w^) - V(A sin2 9)] = 0, 



(r-l)(^-m) 



n + {n — m){2 
and relations 



|V$|2 



-1 



1 



(15) 



Equation (|14| ) with definition (|l^) and relations (|g) and (13) contains actually only flux function ^{r,6), two 
constants s and E as well as metric coefficients. It is elliptic in the subsonic region D > and hyperbolic in the 
supersonic region D < 0. For the details of derivation of this equation see the work by Beskin & Pariev (1993). 
We need to find smooth solution of the equations (p^, (p^, (||) and ( p^ for the function <i>(r, 6*), considering 
boundary conditions at infinity. 
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3 Solution in the supersonic region 

Let us first find tlie solution in the region far under the sonic surface r = for r ^ r^, where we will treat 
it as being close to the solution for radially infalling dust. In the case of cold matter k{s) is equal to 0, so 
II — E = m. By direct substitution into equation (|lj), expressing n in terms of $ by means of the relation (p^), 
it is possible to convince yourself that any arbitrary function $ = $(0) is the solution. This just means that the 
dust falling down to the black hole follows the trajectory 6 = const. At first, we need to specify this arbitrary 
function in order to determine zero order approximation for In our case of nonrelativistic temperature of the 
gas at infinity, sonic surface, at which the local sound velocity is equal to the velocity of the flow, is situated 
far from the black hole horizon th = -M + ^/Ai^ — a? , in the region where relativistic effects are weak. Hence, 
the picture of the stream lines in the transsonic region only slightly deviates from that for the accretion onto 
Newtonian gravitating centre, and the main term in the expansion of $ corresponds to spherically symmetrical 
flow, that is proportional to 1 — cos 61. We will show in section 4 that the corrections to spherically symmetrical 
part of $ in transsonic region are of order of the {M. /ts ) , while the corrections to angular dependent part are 
proportional to [M/ts)'^ and {a/rsY- 

We define small parameter of the expansion e as the difference between the entalpy per particle at infinity 
and the rest particle mass m 

Moo = + e)- (16) 

Far from the sound surface and the black hole u = 0, hence, E = ji. This gives us exact expression for the 
constant E 

E = m{l + e). (17) 
Introducing particle density at infinity Uoo, the expression for entalpy (^ is rewritten in the form 

n 



Tin 



em — . (18) 



It is well known that the accretion rate is determined by the condition of passing the solution for n{r) smoothly 
through the sonic surface (Shapiro & Teukolsky, 1983). Consideration of spherically symmetrical Newtonian 
accretion problem (e.g. Shapiro & Teukolsky, 1983) results in the following accretion rate using our notations 
for entalpy (|l|) and (^ 



$00 = -27rnoo7WV(r) (20) 



where 
and 

,n(r\ - , , , 

2) V 4 y (r- 1)3/2 

This expression is valid for 1 < F < 5/3. The case F = 5/3 needs special treatment and we will not concentrate 
on it though it is this case that is closest to real description of infalling matter if it is fully ionized. Thus, we 
have completely determined the function $o(^'), which was remained to be free in the case of perfectly cold 
matter. We accept the expression (|^ as a zero order approximation to the solution. Note that, when e ^ 
and rioo = const, the accretion rate will become infinitely large. Further, by substitution the expression ([l9|) 
into equation ( p^ zero order approximation for particle concentration no can be found 

no{r) = n^^^{T)e~^/^ (22) 

Also in zero order fiQ — m, Dq = — 1. Radius of the sonic surface is 

.-„.. = Mii^, (23) 



We see that tqs ^ M for small e indeed. 
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The main equation of the problem ( |l4| ) may be expanded in parameter e only when all values, characterizing 
the flow, deviate slightly from those for dust accretion. Therefore, matter needs having nonrelativistic temper- 
ature not only at infinity but also in the whole space. Let us check it. Substitution in expression ( p^ for n its 
zero order approximation (|2^) allows to find entalpy with the first order accuracy 

(r-i)/2 

(24) 



^-m = m(^(r))^-^(5-3r)/2 



2r(r2 



We see that if F < 5/3, for e and r ^ A4 the correction to /i is small indeed and the expansion in e is 

3 /"s" 

now becoming full substantiated. For F = 5/3 radius of the sonic surface ros = "^4^2" ^^^^^ niuch larger 
than A4 , but the temperature in the vicinity of the event horizon is relativistic (for Schwarzschild black hole 
Th ~ 0.111 m, Shapiro & Teukolsky (1983)) and the solution for dust accretion cannot be used as a zero order 
approximation. 

Let us write for the fiux function and particle density $ = $0 + and n = no + ni, where $1 ^ $0 and 
Til <^ no, and perform the linearization of the equation ( [T^ ) with definitio n (p^ and relation (|l3|), allowing to 
express ni via $1. Finally, after eliminating ni from linearized equation (n4) and doing some algebra, stream 
equation can be reduced to the following simple equation in $1 



d^^i 9$i 3r2 



Qj.2 Qj, 2r{r^ + a^) 



-$ooe^"^^^^cos6lsin2 6' 



r-i 



2r(r2 + a^) 



(r-i)/2 



(25) 



The remarkable fact is the cancellation of the derivatives 



and 



o/)9 ' 0/1 o 0/1 ' which makes equation 

to be an ordinary differential equation and enables to find all its solutions. Note, that the equation (Ea) does 
not contain factor A, which vanishes at the event horizon th- This factor is cancelled out from all addenda of 
linearized equation (0), hence, the equation (^5|) has no singularities at the event horizon. There is no even 
the characteristic scale in the latter equation, though from physical point of view one can naturally expect 
it to be revealed. The only characteristic size is a, which can be much smaller than Ai. 

We are able to find easily general solution of equation ( p5| ) , which can be written in the form such that the 
angular dependent part of $1 tends to for r ^ 00 



'I'l = $00 a cos 
+00 



0sin2 2-(r-i)/2(^^P))r-i_^(3r-5)/2^i-3r/2 ^ 



/ 



dr' 



r -\-o 

I [r"(r"2 + a2)]r/2 " ^^^^^ / 



dr' 



i y^r'(r'2 + a2) 



(26) 



where Ai and A2 are arbitrary functions in 9 only. In order to match this solution with the solution in the 
transsonic region r ^ obtained in section 4 let us expand expression (p6|) in a/r ^ 1 



$1=2 



Ai{0) 



1 - 



10 7-2 



24 r4 



$ooCo(-£^) cos0sin2 0((^(F)) 



3r/2 23/ 



2-2r 



M 

_^y3r-i)/2 



00 



]\os9 sin2 9 (^(F))i^-i2-(i^-i)/2ei-3r/2 



(3F-2)(3F- 1) 



(27) 



1 



3F - 1 a2 3 3F - 1 a'^ 



3 + 3F r2 8 7 + 3F r4 

where constant Co denotes the value of the integral 

+00 

Co = 



A2{9), 



47W^ 



-r/2 



d^. 



It is also useful to find the ratio of the perturbation $1 to zero order approximation ^ooe 
leading terms in /r"^ in expansion (p7[) one can obtain 



-3/2 



Retaining only 



1L^3/2 



Co cos^sin^ 9 (^(F))r-i 23/2"2r ,(5-3r)/2 



^lW^3/2 
$00 
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\MJ \ r J ^00 



(28) 



where we have denoted for brevity 



^•"•'-'^"•»^"^"^"'^' (3r-2)(3r-i) ^ (2!.) 

We cannot now determine the functions ^1(6*) and A2{9). It will be done in section 4, where we obtain the 
solution of the stream equation jl^ ) in the region r ^ M/e and match it to (p^ ). At last, it is necessary to 
verify that 4>i/'I>o ^ 1- Indeed, this is immediately seen from the equation (|28D unless functions Ai and A2 
have too high values, which is not the case as it will be shown in section 4. 



4 Solution in the transsonic region 

Now let us consider accretion flow far from the event horizon but still on the scales of order of the radius 
of the sound surface r^. In this region zero approximation for particle density ( p2[ ) can no longer be valid 
because when r ^ JVl/e main term in the difference w^E"^ — fi^ Asrn^ 9 standing for particle density in ( p^ ) 
r^E^ sin^ 9 — jj^r^ sir? 9 vanishes for cold matter, when E' = /i = m, so to determine this difference correctly 
one must take into account thermal corrections to E and /it. Clearly, the consideration of the lowest order 
corrections will lead us to Bondi solution for the accretion in Newtonian gravitational potential of a point mass. 
We have already known the flux function for that solution <I>o given by formula (|9|). Therefore, zero order 
approximations for $ in supersonic and transsonic regions are one and the same and are equal to ^o- However, 
zero order approximations for entalpy and particle density are different. 

To perform expansion of the equations (|l3|), (|l8| ) and ( [T^ ) with definition (^5|) in small parameter e in 
transsonic region, first, we shall introduce some notations. Let us denote the ratio of the particle density to 
that at infinity as y 

y(r^9) = . 

Further it is necessary to introduce rescaled radius x. We define it as 

r = M-^ (30) 

Then, we will search for series in e for y and $ 

y = yo{x) + yi{x, 9) + y2{x, 9) + 7/3(2;, 9) + ..., (31) 

$(r, 9) = $ooe"^/'(l - cos0 + h{x, 9) + ^(a;, 9) + h{x, 9) + .. .), (32) 

where indexes 0, 1,2,3,... mark successive orders of expansion. For n and $ dependent on 9 sonic surface where 
D — will no longer be spherical . Similar to the expansions of y and $ we write an expansion of the radius of 
the sound surface rs{9) 

Xsi9) = xos + xu{9) + X2s{9) +xzs{9) + ..., (33) 



where the zero order value corresponds to ros given by the expression (^3|) 

Because an asymptotical expansion of metric coefficients giu for r — s- 00 contains all successive orders in e it 
is naturally to expect that yi_ and /i will be of order of e, 2/2 and /2 of order of e^, 2/3 and of order of , 
and so on. In further calculation dealing with series we shall follow this order of values that finally will have 
occurred to be true. The only characteristic size in the transsonic region is ros, therefore, 9- and x- derivatives 
of the smooth functions yo, yi, 2/2, . . . , and /i, /2, ... are of the same order and we shall regard in expansions 

~ ~ and so on. In all subsequent expressions indexes 0, 1, 2, 3, . . . mark successive orders of 

ox _a9 ox o9 _ ^ 

expansion of the values considered, index s means that the quantity is taken at the nonperturbed sonic surface 
x = xos, particularly, jj, = fiQ + jii + 112 + ^^^ + ■ ■ D = Dq + Di + D2 + . . ., and similarly for other values. 
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Now zero order approximation for entalpy is obtained by substituting yoi^x) for n/rioc in the expression ( p^ ) 

IIq — m + emyQ~^ . (35) 



1 r-i -L , f ^ /ORN 
1=% -- + :^72Z4' (36) 



Expansion in e of the relation (|13|) to the first order using (|35|), (17), and ( |20D results in the following expression 

r-i 1 I ^'(r) 

which determines the function yo{x). Differentiating ( p^ ) one can obtain 

da; X {T-l)y^+^x^~Lp^Ty ^ ' 

At the point x = xqs, where the denominator in r.h.s. of (|37| ) vanishes, the solution yo{x) must be regular. 
This means that at x — xq^ the numerator must vanish too. Hence, we have three relations among values Xqs, 
yos = yo{xQs), and (p(r); they are equation (|36|), the equality to zero of the denominator in (37) 

<^2(r) = (r-i)%r+i4„ (38) 

and the equality to zero of the numerator in ( ]3^ ) 

<^'(r) = l^lyl. (39) 



Solution of these three relations gives the value of (p{T) just equal to that given by formula (|2l|) used in section 3, 
the value of xqs equal to ( ^^ , and the following value for j/os 

i/(r-i) 



yos = ( ) . (40) 



With the relations (|3^), (|39|) and ( [40| ) the equation ( |36| ) determines two functions yo(a;) passing smoothly 
through the sonic point a; = xqs, one corresponding to Newtonian accretion onto gravitating point object, the 
other to Newtonian ejection from gravitating object (Shapiro & Teukolsky, 1983). Asymptotic of the accretion 

solution for a; — > cxd is ?/o ^ 1 H 1 of the ejection solution j/q ^ x"^; asymptotics for a; — > are: for 

r — 1 a; 

the accretion solution yg — ^x^^^^, for the ejection solution j/o ^ x^^/''^^^\ We plotted functions yo^x) 
on Fig. 1 for T = 4/3. As it can be seen from the asymptotical behaviour of these two functions their relative 
location for all values 1 < F < 5/3 is the same as on Fig. 1. In order to determine the derivative — ^ we 

expand the relation ( |3^ ) up to the second order in the vicinity of the point x — xqs, y — yos- Coefficients in the 
first order expansion identically vanish due to the conditions (38) and (^9|). As a result, we obtain quadratic 
equation in the derivative sought, which has two solutions, one corresponding to the accretion curve yo{x), the 
other to the ejection curve of Bondi solution. Note that no additional conditions should be imposed in order 
to ensure the existence and continuity of all higher derivatives of yoix) at the sonic point. Because r.h.s. of 
the equation ( ^6|) is function analytic in both x and yo in the neighbourhood of a; = xqs two functions yoix) 
determined implicitly by the equation ( ^6|) are analytic too. Since the lines yn — yn (x) are lines of constant value 
of the r.h.s. of the equation (^6|), they cannot intersect each other save the sonic point. Therefore, deducing 
from the asymptotical behaviour of the accretion and ejection solutions, one must choose for the derivative of 
yo{x) at the sonic point larger value from two. Thus, we obtain 



dyo 
dx 



r - 1 / o \r/(r-i) 



x=xas 

Zero order of D given by the formula (|l3) is 



= ) (-4 + yT036r). (41) 



Do = -l + ^/o^'-'- (42) 



At the sonic point Dq{xqs) = 0. In the stream equation (|Tj) 

and substitution of zero order values in all other terms of equation ( p^ ) leads after cancelling out the factor cos 
to the relation (p6|). Hence, in the zero order approximation stream equation (p^ is satisfied automatically. 
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4.1 First order corrections 

Now let us turn to the first order corrections. From (|l^) we obtain 

^^=emyl-^y^{V~l). (43) 
2/0 

Expansion of the quantity entering both the relation ( [l^ ) and the equation ( [T^ ) up to the first order is 



|V$| 



|V$| 



+ (44) 



X2^2(r)sin2 6i V 2/0 sine* 961 



Because of the vanishing of the leading term ~ l/e^ in the relation (jl^) we need to expand it to the second order 
and, hence, one might expect that the rotational parameter a will be involved in the first order approximation, 
because it enters quadratically into the expressions (||) for w'^ and A. However, calculation shows that the 
terms containing are cancelled out from the both sides of the equation ([T^), thus resulting in the following 
expression for yi through /i 

„2/i A„ r-i ^!(r)A , o 2r-2 . V^i^) 2 dh 

2- (^(r- 1)^0 --^j =.(1 + 3^0 -Ay, )--^__ (45) 



Further, we have to linearize equation (|TJ) using (^3|) and ( p4| ) and insert into it the expression ( p5| ) for j/i. The 

^2 

6'-component of the expression embraced by the square brackets in the last term of l.h.s. of (fl3) r 

do 

2Asin0cos0 must be expanded to the second order (because of the vanishing leading term in this difference) 
and, again, a? is cancelled, and the linearization of the stream equation does not contain a?. For D it is 
enough to use its zero order Dp; in the linearized last term in l.h.s. of (|lj) yi stands only multiplied by the 

factor — [(r — 1)?;q^^ — J ^ \ that allows it to be directly replaced by the r.h.s. of expression (p5|). This 
2/0 \ X Ho J 

replacement results in the cancellation of the free term in stream equation. Finally, dividing (|lj) by $oo£ '^^^ 
we obtain the following linear second order differential equation in /i 

Do.'^4^ + {Do + 1) sin4 ( -^M) - ("2. - ^) ?A . o, (46) 



where is given by the expression (p2[). 

Equation (Em is an equation with separable variables and angular operator 



L,=sin0|-('-i-|-) (47) 



de \smede 

is the same as in the linearizations of the stream equation considered in BP and by Petrich, Shapiro & Teukol- 
sky (1988). Particle flux density must be finite at the rotational axis. As it is seen from the expression (0) for 

^ dfi 

?iup this leads to that the quantity — must be finite at 6* = and 9 ^ tt. Besides, it is obvious that at 

sin ay 

^ = total flux must be equal to 0, so one more boundary condition for /i is fi\e=Q — 0. It can be proved that 
the eigenvalue problem for operator Lg with these boundary conditions has eigenvalues Cn — —n(n+ 1), where 
n = 0, 1, 2, 3, . . ., with corresponding eigenfunctions 

Qn= „ '^- sin^gPJcosg), 

(2n)! 

where prime denotes differentiation of the Legendre polynomial PniO with respect to ^ and normalizing factor 
is chosen the same as in BP to make the comparison of our results to the results by BP more direct. First four 
eigenfunctions are 

Qo^l-cos^, Qi=sin2 6', Q2 = sin^ 6* cos 6*, Q3 = sin^ 6* f cos^ 6* - ^ j . (48) 
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We look for the solution of the equation ( p6| ) in the form 

oc 

fl^Y.^9ln{x)Qn{0). (49) 



ji=0 



Then gin{x) must satisfy the following equation 



Dox^^ + ( 44, ~2x)^- {Do + l)n(n+ l)5i„ = (50) 



dx"^ \(p^(T) J dx 

with boundary conditions that gin{x) must be regular at the sonic point where -Do = and be finite at a; — s- oo. 
Two linearly independent solutions of ( |50| ) have asymptotics ^ a;"+^ and ~ so we must choose the latter 
solution. Because of the relation ( p9| ) it follows from equation ( |50| ) that n(n + l)gin{xQs) = 0. For n 7^ it 
means that ginix) = 0, while for n = one obvious solution of uniform equation ( p0| ) is gw — const. Second 
linearly independent solution has asymptotic gio x for x 00 and must be ruled out. Thus, first order 
correction to the flux function fi has the form 

A = e<?io(l-cos0). (51) 

To complete construction of the first order approximation it remains to calculate j/i according to the equa- 
tion ( ^ . We see that yi does not depend on 9. Correction to particle density yi {x) must be smooth function in 
the vicinity of the sonic point x = xqs, where according to ( ^ l.h.s. of the equation (^) vanishes. Therefore, 
r.h.s. must vanish too, which allows us to determine the value of gw 

510 = 7 K (52) 



Now j/i is given by the formula 



4 5-3r 



1 I Q 2r-2 A r-i o (r) 
1 + 3% -4yo -2510^^^ 

= eyo r::^ — r-rri — — -■ (53) 



2(r - l)yr-i _ 2 



Both numerator and denominator in ( |53| ) are analytic functions in yo and x, therefore, yi{x) is also analytic 
function in the neighbourhood of the sonic point, and no additional constrains should be imposed in order of 
continuity of the derivatives of yi . The value of yi at the sonic point is 



yi \x=xo. 



2(5-3r)vio-6r 



2ir- 5 
r + 1 



-4 + Vio-er + 6(3r + 1 



(54) 



Asymptotical behaviour for x cx3 is j/i ^ — -. Function yi(x) is plotted on Fig. 2 for some values of F. 

. ^(^~^'' 

What we actually have obtained considering first order approximation are the corrections to Newtonian 
accretion rate and particle density in the problem of spherically symmetrical accretion onto a Schwarzschild 
black hole. As long as a is not contained in the results they coincide with those for a Schwarzschild black hole. 
Indeed, exact analytical solution for the accretion onto a Schwarzschild black hole (e.g., Shapiro & Teukolsky, 
1983) in our notations can be written as 

$ = <^>sw{l - COS0) 



with 



(r-l + ^(2-3F) 



i/(r-i) 



(55) 



M 

where the value ~ e should be obtained by solving cubic equation 

2r. 
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For the rescaled radius of the sonic surface Xgw = tTsj M. expansion of the sohition of equation ( |56| ) in terms 
of e up to e' is 



5-3r 

4(r-i) 



3or-9r2-i3 2 (sr^ ~ i4r + 3)(r - 1) 



2(5-3r)2 



(5 - 3r)'^ 



e^TTT^ — rT^(3r + i)(9r3 - Tsr^ + losr - 53) 



8(5-3r)6 



(57) 



For all 1 < r < 5/3 the successive terms in the series ( p7[ ) have alternating signs: Xsw\ > 0, Xsw2 < 0, Xsw3 > 0. 
Substitution of the expression (57) into formula (|55|) results in series for the accretion rate 



1 



3 3r + 1 



3 135r3 - 99r2 - 207r + 139 



4 5 - 3r 32 (5 - 3r)3 

3 8505r5 - 9153r4 - 37530r3 + 89226r2 - 71583r + 20279A 



128(5- 3r)^^ 



(58) 



First order term in (|5^) coincides with the expression for 510 (^2|) and is positive, which means enhanced 
accretion rate. Term proportional to ig positive for 1 < F < 1.3037 and negative for 1.3037 < F < 5/3, term 
proportional to is positive for 1.1558 < F < 5/3 and negative for 1 < F < 1.1558. 

At last, let us show that the correction to the radius of the sonic surface xis (see expansion (p3|)) does not 
depend on and is equal to Xswi from (|57|). In general sonic surface is located where D = 0. First order of D 
given by the formula (15) is 



Di = iT~ 1) 



(p2(r) 



yo X 



(59) 



and does not depend on 9. Correction xis is determined from the equation 

dDo 



dx 



Di{xos 



0. 



Its solution is readily obtained using expressions ( |4l| ) and ( p4| ) 

30F - 9F2 - 13 



Xls 



8(F-1)(5-3F) 



(60) 



and coincides with Xs 



4.2 Second order corrections 

In order to match the solution in the supersonic region we need to consider further orders of expansion in the 
transsonic region. Finding 1/2 and /2 we will act similarly as when we deal with yi and /i. Let us denote by 
d — a/Ai dimcnsionless value of a. From (|l8|) we obtain 



r-i/-p T \ I 

m = emy^ (F - 1) — 



r-2yl 



vl 



further, 



|V$|' 



7^2(^2 (r) sin2 



2 

Vl 



■ cos2 Oyl - iegioyoyi 



2 I o 2 2 2 



de 



Third order of expansion of relation (|13D gives an expression for ?/2 as a function of unknown /2 



,y2 



VQ 

-510(510 



(F-l) 



x'v'o 



ra2 <^2(p) 



a;2 ygX'^ 



2— cos^ 
x^ 



vl-') 



Vq- 



2e 



yo 



yl^-\T-l)-2{T-l)y^,-' 



^vf-'ivl-' 

Vlx' 



+2g 



Vlx^ 



vl 
vl 



(F-l)(F-2)2;i 



r-i 



ylx^ 



2 dh^\v) 



sin 61 dQ 



Vo^^ 



(61) 



12 



Note that when obtaining second order of the stream equation (|T^) we need to retain only leading term Dq 
in the series for D because both (f>o and $i do not depend on r and Le$o = Lg^i = 0. Therefore, partial 
differential equation for /2 changes its type at the surface x = xqs as well as equation ( p6[ ) for /i does. Again, 

y2 appears in the second order of the stream equation only multiplied by — T j/Q^^(r — 1) — , which 

allows it to be excluded by means of expression (^l]). Finally, dividing by <i>ooe^'^^^ one can obtain the 
following linear second order partial differential equation in /2 

Dox ^ + {Do + l)iiine— {-——]- [2x - 



e2^cos0sin20(l-^^|^) =0. (62) 



dx^ ^ ' de Vsin6i dO J \ f'^iV) J dx 

■^cos6'sm^6' 1 ^—r 

x'^ \ xi-p'^yy)^ 

Variables in equation (^2|) can be separated and angular operator coincides with that in equation (p6|). Moreover, 
equation (p2h differs from equation (^6|) only by that the former contains free term, while the latter does not. 
Equation (p2) coincides with the limit ^ 1 of the equation (52) in BP, which is the linearization of the 
stream equation in small parameter aj M.^ so the results of BP concerning the case <g; 1 can be applied to 
the equation (^2|). Similar to the equation ( ^ ) we seek for the solution in the form 

oo 

Then, radial functions g2n{x) for all n 7^ 2 must satisfy exactly the same equation as ( |50| ) with identical boundary 
conditions. This immediately allows to conclude that g2n = for all rt 7^ 0, 2 and 1720 = const. As for 322 it 
must satisfy the following equation 

which is identical to equation (57) in BP. The requirement of regularity of 522(2;) at nonperturbed sound surface 
Dq — leads to the condition 



g22{xos) = -TTv-- (64) 

Another boundary condition for 522 is 522 ^ x^"^ when x ^ 00. Similarly to BP we introduce function G{x) 
such that the solution of (E3) satisfying boundary conditions is written in the form 



522(2;) = -G{x)- 



^ , (i-3r)/2 

XOs 



where G{xos) = 1/2 and G{x) - a;(3r-5)/2 fo^. 

X — > 00. The dependencies G{x) found numerically by integrat- 
ing ( |63| ) with boundary conditions are shown on Fig. 3, which is actually the same as Fig. 1 in BP. Thus, second 
order correction to the flux function /2 has the form 

-2 / \ (i-3r)/2 

/2 = e^52o(l-cos6i)-e2G(a;)-^ I — ) sin^ 61 cos 61. (65) 

Xqs \XQsJ 

Similar to constant 510 constant 520 should be determined from the requirement of smoothness of the value yi 
at the sonic surface x = xqs- Inserting into expression ( |6l| ) value for /2 from (^) we see that 

^ ^-^^ e2g2O + e'.922(a;)(3cos2 0- 1). 



sin 6* de 



Therefore, expression ( |6l| ) contains terms proportional to cos^ 6 and terms dependent only on x. For evaluating 
r.h.s. of (^ at the nonperturbed sonic surface we use expressions (64) for 522(a;os), (jsj) for yi{xos), (^ ) 
for yo{xos), and ( ^2[ ) for 510. As a result, terms proportional to cos^ vanish, while the requirement for the rest 
addenda to be equal to determines the value of 520 
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which actually coincides with the second order term in the expansion of ([S^) ■ This result is obvious because 
in the case = when 522 = and /2 = £^.920(1 — cos 9) we deal with the accretion onto a Schwarzschild black 
hole. Now the construction of the second order approximation is completed. 

Now we are able to match the solution in the supersonic region obtained in section 3 with the solution in 
the transsonic region. For this purpose we should consider the behaviour of 1722(2;) for a; — > 0. Expanding the 
relation (36) in the vicinity of x = one can find two leading terms in the asymptotical series for yo{x) 



yo 



V2- 



-3/2 



1 



r-i 



2(r+i)/2 



.(5-3r)/2 ^ p |^^(5-3r)/2^ 



Substitution of this expression into r.h.s. of (|63[ ) allows us to find first nonzero term in the expansion of r.h.s. 
of (|6^). Hence, in the limit x ^ equation (|63|) takes the form 



(fg22 , 3 dg: 



22 .2 

= —a 



(^(r)) 



r-i 



dx"^ 2x dx 
General solution of the last equation is 

522 = Cix-^/^ + C2 



2(r-i)/2 



„-3(r+i)/2 



^2^(l-3r)/2 



(67) 



where the constant K{V) occurs to be just equal to that given by (|29|) and constants Ci and C2 are uniquely 
determined by boundary conditions for (722- Since 522(2:) satisfies boundary condition (|^ at the sonic point, 
constants Ci and C2 are of order of and it is convenient to write them explicitly Ci — ci(r)a^, C2 — C2iT)a? . 
Therefore, in the inner region x ^ e the third term in the expression ( |6^ ) is dominant. This conclusion is 
confirmed by numerical computation of (722(2;), the results of which are plotted on Fig. 3. At a; ^ function G{x) 
tends to a constant value, which demonstrates that the first two terms in the expression (67) become negligible 
compared to the third one. Limiting value of G[x) is derived from the last term of ( |67| ) 

hm G{x) = 



(3r-2)(3r-l)(r-l)2r 



and is in good agreement with the numerical results. Limit r S> of the inner solution ( |28|) when expressed 
in terms of rescaled variable x takes the form 



1l 



,3/2 



$oov 



a cos a sm 



For outer solution (32) combining ^ 



= 2CqS? cos e sin2 0(^(r))i^-i23/2-2r2,-i/2,3(2-r)/2 „ 2 ^^L a^'^/^e^ 



2 0i^(r)x(l-3r)/2^2^^^3/2_ 



we have 



(68) 



1l 



,3/2 



e^Ci cos e sin^ Bx-^'^ + e^C2 cos sin^ - e^d'^x'^^-^^^^^K{T) cos sin^ 



+ (1 - cos6')(e5io + 6^520) + 0{e^), 



(69) 



where symbol 0(e") denotes terms proportional to e". In order to match expression ( pSj ) to (|69| ) we must choose 
the following values for functions Ai (0) and A2 (6) 



Ai (0) = <S>aoVMa^ cos sin^ 9 



Co(^(r))i^-i2^/2-2r,i-3r/2_ l^^(r) 



Oie), 

A2{0) = $ooa^e^^^C2(r) cos 61 sin2 + $00(1 - cos 6l)e-i/2(.gio + e.g2o) + 0(e^/^). 

With these expressions solution (|2^) in the supersonic region becomes fully determined and can be expressed 
by the following formula 



-\-oo 

^£3/2 ^ _^cos0sin2 0((p(r))r-i2i-2r,(5-3r)/2 f 



r/{2M) 



v/r(t2+«V4) 



+00 



dc 



a2/4)]r/ 



9 + cos sm 



,ci(r) 



3/2 



r/{2M) 



(70) 



e'-^C2(r)cos(?sin2 + {1 ~ cos9){egio + e'520) + 0{e^/^) 
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Firstly, we see that $1 ^ $o- Secondly, dominating term in (jT^) is the first one and, since we deal only 
with the first order approximation to we cannot be sure that higher order terms are negligible in comparison 
with the rest addenda in ([70|). At the same time, only terms dependent on 9 like Qo ^ I — cos 9 contribute 
to the total accretion rate, and because the total accretion rate must be some constant value independent of 
the radius these terms do not depend on r in all orders in e. Hence, these terms are the same in the supersonic 
region and in the transsonic region and are determined correctly in the formula ( |70| ) up to the 0{e^). Finally, 
we write down explicitly the solution in the transsonic region which matches the solution ( fzol ) 

$ ^ $ooe"''/'[l + (egio + e'.g2o)(l - cos 9) + e^g22{r) cos^sin^ 9 + 0{e^)]. (71) 

On Fig. 4 we plotted the picture of the stream lines corresponding to the flux function given by formula (|70| ) 
for r ^ and by formula ( [fl]) for r ^ r^g. Stream lines concentrate toward the equatorial plane of a black 
hole. 

To find second order correction X2siS) to the shape of the sonic surface one should expand the condition 
D{xs) = up to the second order in e similarly to that was done when determining first order correction xis (|60|). 
However, because the second order solution in the transsonic region is actually the same as the solution for 
a ^ M and Vs ^ M obtained in BP, it is much easier to notice that the shape of the sonic surface up to the 
second order in e coincides with that found in BP for the case a <C and Ai jrg <C 1. So, rewriting formula (74) 
from BP and adding term proportional to from expression for Xsw (|57|), we obtain in our notations 



e 



7W2(5-3r)2L^ ' y J jy /v ;j (5 - SF) 



where fc,(F) ^ 4s^ ^ 



is computed numerically in BP (see Table 1 there), ^2(4/3) = 0.861. The 



X — XQs 



shape of the sonic surface is also shown on Fig. 4 for F = 4/3 and a ~ A4. 



5 Conclusions 

Using stream (Grad-Shafranov) equation (|lj) we have considered the problem of hydrodynamic accretion onto 
a resting Kerr black hole. The fact that in astrophysical situations, when this type of accretion could be 
realized (accretion of interstellar medium onto an isolated black hole), the temperature of the surrounding gas 
is nonrelativistic introduces small parameter in the problem and allows to obtain analytical solution of the 
stream equation (|l|) in the case of an arbitrary value of the rotational parameter a, extending thus the results 
of the work BP onto rapid rotation of the black hole. 

Form of the sonic surface and corrections to the particle density are determined. However, because the sonic 
surface is located far from the horizon of a black hole these corrections are very small and flow picture does not 
differ drastically from the case of the accretion onto a Schwarzschild black hole save the dragging into rotation 
around Kerr black hole with Lense-Thirring angular velocity. Therefore, the related astrophysics (emergent 
radiation) does not change significantly from previous consideration by Shapiro, 1974. Presence of a magnetic 
field or small angular momentum of the gas at infinity relative to the black hole have much more influence on 
the accretion pattern and emergent radiation. So, the results of the present work lie primarily in the scope 
of mathematical physics and show that shock-free smooth solution exists in the problem of the steady state 
hydrodynamic accretion onto a rotating black hole with arbitrary a < M. 
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Figure Captions 

Figure 1. Functions yo{x) corresponding to accretion and ejection solutions of equation ( p6| ) for F — 4/3. 
Figure 2. Dependencies of yi on x for different values of F. 
Figure 3. Function G{x) for different values of F. 

Figure 4. Stream lines and sonic surface for e — 0.3, a — A4 and F 4/3. 



